Instability and new phases of higher-dimensional rotating black holes 
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It has been conjectured that higher-dimensional rotating black holes become unstable at a suffi- 
ciently large value of the rotation, and that new black holes with pinched horizons appear at the 
threshold of the instability. We search numerically, and find, the stationary axisymmetric pertur- 
bations of Myers-Perry black holes with a single spin that mark the onset of the instability and 
the appearance of the new black hole phases. We also find new ultraspinning Gregory-Laflammc 
instabilities of rotating black strings and branes. 



Black holes are the most basic and fascinating objects 
in General Relativity and the study of their properties 
is essential for a better understanding of the dynamics 
of spacetime at its most extreme. In higher-dimensional 
spacetimes a vast landscape of novel black holes has be- 
gun to be uncovered [T|. Its layout — i.e., the connec- 
tions between different classes of black holes in the space 
of solutions — hinges crucially on the analysis of their 
classical stability: most novel black hole phases are con- 
jectured to branch-off at the threshold of an instability 
of a known phase. Showing how this happens is an out- 
standing open problem that we address in this paper. 

The best known class of higher-dimensional black 
holes, discovered by Myers and Perry (MP) in [2], ap- 
pear in many respects as natural generalizations of the 
Kerr solution. In particular, their horizon is topologi- 
cally spherical. However, the actual shape of the horizon 
can differ markedly from the four-dimensional one, which 
is always approximately round with a radius parametri- 
cally ~ GM. This is not so in d > 6. Considering for 
simplicity the case where only one spin J is turned on (of 
the |_^tM independent angular momenta available), it is 
possible to have black holes with arbitrarily large J for a 
given mass M. The horizon of these ultraspinning black 
holes spreads along the rotation plane out to a radius 
a ~ J/M much larger than the thickness transverse to 
this plane, r+ - (GM 3 / J 2 ) 1 /^) _ xhis fact was p i c k e d 
out in [3] as an indication of an instability and a con- 
nection to novel black hole phases. In more detail, in the 
limit a->oo with r + fixed, the geometry of the black hole 
in the region close to the rotation axis approaches that 
of a black membrane. Black branes are known to exhibit 



classical instabilities [4 , at whose threshold a new branch 
of black branes with inhomogeneous horizons appears [5] . 
Ref. [5] conjectured that this same phenomenon should 
be present for MP black holes at finite but sufficiently 
large rotation: they should become unstable beyond a 
critical value of a/r + , and the marginally stable solu- 
tion should admit a stationary, axisymmetric perturba- 
tion signalling a new branch of black holes pinched along 
the rotation axis. Simple estimates suggested that in fact 
(a/r + ) cr it should not be much larger than one. As a/r + 
increases, the MP solutions should admit a sequence of 
stationary perturbations, with pinches at finite latitude, 
giving rise to an infinite sequence of branches of 'pinched 
black holes' (see fig. [IJ. Ref. [B] argued that this struc- 
ture is indeed required in order to establish connections 
between MP black holes and the black ring and black Sat- 
urn solutions more recently discovered. Our main result 
is a numerical analysis that proves correct the conjecture 
illustrated in fig. [I] 

The solution for a MP black hole rotating in a single 
plane in d dimensions is [2] 

r d-3 /J r 2 \ 

ds 2 = -dt 2 + -J^(dt + a sin 2 Od<p)~ + Z( — + d0 2 \ 
+(r 2 + a 2 ) sin 2 6 d<p 2 + r 2 cos 2 dft 2 d _ 4) , (1) 

£ = r 2 + a 2 cos 2 fl, A = r 2 + a 2 - ^ . (2) 

The parameters here are the mass-radius r m and the 
rotation-radius a, 

d _ 3 _ IGttGM _ d — 2 J 

Tm ~ (d-2)Sl d - 2 ' a ~ ~ 2~ M ' ( ) 
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FIG. 1: Diagram of entropy vs. spin, at fixed mass, for MP 
black holes in d > 6 illustrating the conjecture of ;3] (see also 
[6]): at sufficiently large spin the MP solution becomes un- 
stable, and at the threshold of the instability a new branch 
of black holes with a central pinch appear (A). As the spin 
grows new of branches of black holes with further axisymmet- 
ric pinches (B, C, . . . ) appear. We determine the points where 
the new branches appear, but it is not yet known in which di- 
rections they run. We also indicate that at the inflection point 
(0), where d 2 S/dJ 2 = 0, there is a stationary perturbation 
that should not correspond to an instability nor a new branch 
but rather to a zero-mode that moves the solution along the 
curve of MP black holes. 



The event horizon lies at the largest real root r — r + of 
A. 

The linearized perturbation theory of the Kerr black 
hole (d = 4) was disentangled in [7] using the Newman- 
Penrose formalism. Attempts to extend this formalism to 
decouple a master equation for the gravitational pertur- 
bations of (JT|) in d > 5 have failed so far. Moreover, even 
if some subsectors of the perturbations of some classes 
of MP black holes have been decoupled [8] , none of them 
shows signs of any instability and indeed they do not con- 
tain the precise kind of perturbations we are interested 
in. Thus we take a more frontal numerical approach to a 
full set of coupled partial differential equations (PDE) . 

We intend to solve for a stationary linearized perturba- 
tion h ab around the background . Choosing traceless- 
transverse (TT) gauge, h a a = and V a h ab = 0, the 
equations to solve are 

(A L h) ab = -V c V c h ab - 2R a c b d h cd = , (4) 

where is the Lichnerowicz operator in the TT gauge. 
Actually, we solve the more general eigenvalue problem 

(A L h) ab = -k 2 h ab , (5) 

which is known to appear in two contexts: eqs. |5]) de- 
termine the stationary perturbations of a black string in 
d+1 dimensions (obtained by adding a flat direction z to 
([lj) with a profile e lkz h ab . Thus such modes with k > 
correspond to the threshold of the Gregory-Laflammc in- 
stability of black strings [I]. The same equations also 
describe the negative modes of quadratic quantum cor- 
rections to the gravitional Euclidean partition function 



[S]. A recent study of this problem for the Kerr black 
hole has shown the existence of a branch of solutions ex- 
tending the negative Schwarzschild mode (with fcg c h 7^ 0) 
to finite rotation, with k growing as the rotation increases 
towards the Kerr bound [TP] , 

Our reason to consider ([5| instead of trying to solve 
directly for k = is that there exist powerful numerical 
methods for eigenvalue problems that give the eigenval- 
ues k together with the eigenvectors, i.e., the metric per- 
turbations. If the ultraspinning instability is present for 
MP black holes in d > 6, then, in addition to the analogue 
of the branch studied in [10] . a new branch of negative 
modes extending to k = must appear. The eigenvalue 
k = corresponds to a (perturbative) stationary solution 
with a slightly deformed horizon. In fact, as explained 
above, we expect an infinite sequence of such branches 
that reach k = at increasing values of the rotation. The 
solutions for k > imply new kinds of Gregory-Laflamme 
instabilities and inhomogeneous phases of ultraspinning 
black strings (see also [TT]). 

The modes we seek preserve the SO(d — 3) x SO (2) 
rotational symmetries of the MP solution and depend 
only on the radial and polar coordinates, r and 9 [3J. 
Thus we take the ansatz 

ds 2 = -e 2 " (dt-cud(f>) 2 + e 2,/1 d<p 2 + e 2 ' n sin 2 9 dd 2 

+e 27 (dr - x sin 9d9) 2 + e 2 *dO^_ 4 . (6) 

We decompose a given quantity Q — {vq, v\, 77, 7, w, \} 
as Q = Q + SQ. The unperturbed contribution Q{r, 9) 
describes (JT|). The perturbations SQ(r, 9) are determined 
solving the eigenvalue problem (|5| subject to appropriate 
boundary conditions. After imposing TT gauge, eq. §5§ 
reduces to four coupled PDEs for Sij, dj, 8\ and 5$ (the 
TT conditions then give Suq, 8v\ and 80S). The bound- 
ary conditions are that the perturbations are regular and 
finite at the horizon, r = r + , at infinity, r = 00, and at 
the poles 9 — 0, tt/2. In addition, we impose 8x(r+) = 0. 
It is important to ensure that the eigenmodes we find 
are not pure gauge, h ab = V( a ^&). We can prove that 
in the TT gauge, pure gauge perturbations within our 
ansatz necessarily diverge at either the horizon or infin- 
ity. Thus, with our boundary conditions, the eigenmodes 
we obtain are never pure gauge. 

We use a numerical approach successfully applied to 
the identification of the negative mode of Kerr and Kerr- 
AdS black holes [TU]. It employs a Chebyshev spectral 
numerical method (see [TO] for further details). We have 
carried out the calculations for d = 7, 8, 9. The cases d = 
5 (where the heuristics of [3] do not allow to predict any 
instability) and d — 6 present more difficult numerics. 
These, as well as a more detailed presentation of our 
numerical approach, will be discussed elsewhere. 

The results for d—1 are displayed in fig. [2] the other 
two cases being qualitatively very similar. We plot the 
negative eigenvalue — k 2 as a function of the rotation pa- 
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FIG. 2: Negative eigenvalues for the MP black hole in d — 7. 

rameter a. We normalize k and a relative to the mass- 
radius r m , which is equivalent to plotting their values for 
fixed mass (or mass per unit length, in the black string 
interpretation). As described above, the leftmost curve, 
which does not reach k — 0, is the higher-dimensional 
counterpart of the Kerr negative mode, and the eigen- 
values k are the wavenumbers of the Gregory-Lafiamme 
threshold modes at rotation a. At larger rotation we find 
new branches of negative modes that intersect k = at 
finite a/r m . We label these successive branches with an 
integer I = 1,2,3,..., and refer to them as 'harmonics'. 
The values of a/r m at which the stationary perturbations 
appear are listed in table [T] 



d 


(a/r m )\t=i 


(a/r m )\e=2 


(a/r m )\e=3 


7 


1.075 


1.714 


2.141 


8 


1.061 


1.770 


2.275 


9 


1.051 


1.792 


2.337 



TABLE I: Values of the rotation a/r m for the first three har- 
monics of stationary perturbation modes (k = 0). The esti- 
mated numerical error is ±3 x 10" 3 in d = 7 and ±5 x 10~ 3 
in d = 8,9. 

It is important to note that the k — cigcnmode of 
the harmonic i = 1 does not correspond to a new sta- 
tionary solution. Instead it is a zero-mode that takes the 
solution to a nearby one along the family of MP black 
holes. The existence and location of this zero-mode is a 
consequence of the fact that if the Hessian of the Gibbs 
potential G(T, n t ) = M-TS-J2i calculated along 
a family of solutions, has a zero eigenvalue for some par- 
ticular solution, then there is a zero-mode perturbation of 
the gravitational (Euclidean) action / = G/T that takes 
that solution to an infinitcsimally nearby one along that 
family. That is, perturbing the solution with that zero- 
mode does not correspond to branching-off into a new 
family of solutions. 

One can easily check that the determinant of the Hes- 
sian of the Gibbs potential is proportional to the deter- 
minant of the Hessian of the entropy with respect to only 



the angular momenta, i.e., to the determinant of 

Therefore, for solutions with a single spin, there must 
appear a stationary perturbation, in principle not asso- 
ciated to an instability of the black hole, at the inflection 
point of the curve S(J) at fixed M (point in fig. [TJ. 
For the MP solutions this happens at 

The values of (a/r m ) mom for d — 7,8,9 agree with the 
central values of the numerically-determined rotations 
(a/r m ) for t = 1 (first column in table [I]) up to the third 
decimal place. This is a very good check of the accuracy 
of our numerical methods. 

The k — eigenmodes of the higher harmonics, I > 2, 
do not admit this interpretation as perturbations along 
the MP family of solutions and thus correspond to gen- 
uinely new (perturbative) black hole solutions with de- 
formed horizons. Their appearance conforms perfectly 
to the predictions in [5] and [BJ. It is then natural to 
expect, although our approach does not prove it since it 
only captures zero-frequency perturbations, that the har- 
monic I = 2 signals the onset of the instability conjec- 
tured in [3] . The k = eigenmodes for higher harmonics 
confirm the appearance of the sequence of new black hole 
phases as the rotation grows. 

To visualize the effect on the horizon of the perturba- 
tions that give new solutions, and provide further con- 
firmation of our interpretation, we draw an embedding 
diagram of the unperturbed MP horizon and compare it 
with the deformations induced by the ultraspinning har- 
monics I > 2. This is best done using the embedding 
proposed in [12], which has the advantage of allowing to 
embed the horizon along the entire range < 6 < ir/2 for 
any rotation, although at the cost of stretching the pole 
region, which acquires a conical profile. We do it for the 
I = 2,3,4 ultraspinning harmonics in figs. [3] [5] In spite 
of the distortion created by the embedding, the effect of 
the perturbations is clear: I = 2 modes create a pinch 
centered on the rotation axis — 0; I — 3 modes have 
a pinch centered at finite latitude 8; I = 4 modes pinch 
the horizon twice: around the rotation axis and at finite 
latitude. These are the kind of deformations depicted 
in fig. [l] To better identify the number of times that 
the perturbed horizon crosses the unperturbed solution, 
in these figures we also plot the logarithmic difference 
between the two embeddings. 

Ref. [3] gave several arguments to the effect that criti- 
cal values a/r m close to 1 were to be expected. In partic- 
ular, it was pointed out that the change in the behavior of 
the black hole from 'Kerr-like' to 'black membrane-like' 
could be pinpointed to the value of the spin where the 
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FIG. 3: Embedding diagram at (a/r m ) cr lt of the d = 7 black 
hole horizon, unperturbed (solid), and with the first unstable 
harmonic perturbation (i = 2, k = 0) (dashed). The em- 
bedding Cartesian coordinates Z and X lie resp. along the 
rotation axis 9 = and the rotation plane 6 = n/2. We also 
show the logarithmic difference between the embeddings of 
the perturbed (Zt-2) and unperturbed (Zq) horizons. The 
spikes represent the points where the two embeddings inter- 
sect. The perturbation has two nodes, so the horizon squeezes 
around the rotation axis, then bulges out, and squeezes again 
at the equator, as in the conjectured shape A in fig. [T] 
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FIG. 4: Like fig. [3] for £ = 3: between the first two nodes of 
the perturbation the horizon has a pinch (shape B in fig. [I]). 



a = 2.489 r m 




FIG. 5: Like fig. [3] for £ — 4: the four nodes deform the 
horizon into shape C of fig. [T] 



temperature (i.e., surface gravity) reaches a minimum 
for fixed mass, which is the same, for solutions with a 
single spin, as the inflection point of S(J). As we have 
argued, the zero-mode at this solution should not signal 
an instability. The 1 = 2 mode at the threshold of the 
actual instability instead appears at larger rotation, well 
within the regime of membrane-like behavior as conjec- 
tured in [6]. We expect this to be true in general: the 
ultraspinning instability of MP black holes should appear 
for angular momenta strictly beyond the (codimension 1) 
locus in the space of angular momenta where the Hessian 
Hij has a zero eigenvalue. 

In particular, in d = 5 this criterion does not allow 
any ultraspinning instability for any J±, J2, and in d > 6 
with all the N = angular momenta Jj equal it 

predicts that the instability should appear at a/r m > 
2-W/(d-3)_ jj owever we cannot predict the precise values 
of the rotation where the instability appears. 

We have identified the points in the phase diagram 
where the new branches must appear, but we cannot de- 
termine in which direction these run. This requires calcu- 
lating the area, mass and spin of the perturbed solutions. 
However, for any k ^ — and numerically we can never 
obtain an exact zero — the linear perturbations decay ex- 
ponentially in the radial direction, and so the mass and 
spin, measured at asymptotic infinity, are not corrected. 
It seems that in order to obtain the directions of the new 
branches one has to go beyond our level of approximation 
or adopt a different approach. 

The new I > 1 branches extend to non-zero eigenvalues 
k. These imply a new ultraspinning Gregory-Laflamme 
instability for black strings, in which the horizon is de- 
formed not only along the direction of the string, but 
also along the polar direction of the transverse sphere. 
Observe that, even if the £ = 1, k = mode is not 
an instability of the MP black hole, the modes £ = 1, 
k > are expected to correspond to thresholds of GL 
instabilities of MP black strings. At a given rotation, 
modes with larger £ have longer wavelength fc _1 and so 
the branch £ = 1 is expected to dominate the instability. 
The growth of k with a can be understood heuristically, 
since as a grows the horizon becomes thinner in direc- 
tions transverse to the rotation plane and hence it can fit 
into a shorter compact circle. 

To finish, we mention that pinched phases of rotat- 
ing plasma balls, dual to pinched black holes in Scherk- 
Schwarz compactifications of AdS, have been found |13j . 
as well as new kinds of deformations of rotating plasma 
tubes [T3] and rotating plasma ball instabilities [T5] . The 
relation of our results to these and other phenomena of 
rotating fluids will be discussed elsewhere. 
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